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Abstract. We study the statistics of 61 measured masses of neutron stars (NSs) in binary pulsar systems, including 18 double 
NS (DNS) systems, 26 radio pulsars (10 in our Galaxy) with white dwarf (WD) companions, 3 NSs with main-sequence 
companions, 13 NSs in X-ray binaries, and one undetermined system. We derive a mean value of M = 1.46 ± 0.30M Q . When 
the 46 NSs with measured spin periods are divided into two groups at 20 milliseconds, i.e., the millisecond pulsar (MSP) group 
and others, we find that their mass averages are, respectively, M = 1 .57 ± 0.35M o and M = 1 .37 ± 0.23M Q . In the framework of 
the pulsar recycling hypothesis, this suggests that an accretion of approximately ~ 0.2M o is sufficient to spin up a neutron star 
and place it in the millisecond pulsar group. Based on these estimates, an approximate empirical relation between the accreting 
mass (AM) of recycled pulsar and its spin period is proposed as AM = 0.43(M o )(/ > /lmi) _2/3 . If we focus only on the DNS, 
the mass average of all 18 DNSs is 1.32 ± 0.14M G , and the mass averages of the recycled DNSs and the non-recycled NS 
companions are, respectively, 1.38±0.12M G and 1.25+0.13 M G . This is consistent with the hypothesis that the masses of both 
NSs in DNS system have been affected by accretion. The mass average of MSPs is higher than the Chandrasekhar limit 1.44M Q , 
which may imply that most of binary MSPs form via the standard scenario by accretion recycling. If we were to assume that the 
mass of a MSP formed by the accretion induced collapse (AIC) of a white dwarf must be less than 1.35 M Q , then the portion of 
the binary MSPs involved in the AICs would not be higher than 20%, which imposes a constraint on the AIC origin of MSPs. 
With accreting mass from the companion, the nuclear matter composition of MSP may experience a transition from the 'soft' 
equation of state (EOS) to a 'stiff' EOS or even neutron to quark matter. 
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1. Introduction 

Mass is one of the important parameters of a neutron star (NS), 
from which we can infer the stellar evolution of its progeni- 
tor, the nuclear matter composition of a compact object (e.g. 
Haensel et al. 2007) or its equation of state (EOS), and strength 
of gravitational field if the NS radius is known. In other words, 
the precise mass measurements can provide significant tests of 
studies of stellar evolution, nuclear physics of superdense mat- 
ter, and Einstein's general relativity in the strong gravity regime 
(Lattimer & Prakash 2004, 2007; Kramer & Stairs 2008), as 
well as insight into binary evolution since NS masses are mea- 
sured in binary systems. 

A NS is one of the possible ends for a massive star with 
mass greater than ~ 4 - 8 M Q . After having finished burning the 



nuclear fuel, a star undergoes a supernova (SN) explosion, and 
the central region of the star collapses under gravity to form 
a NS in the central supernova remnant (SNR) (Haensel et al. 
2007). Hence, the NS mass statistics help the astronomers to 
infer the properties of its progenitor star, and its links to SN 
and SNR. However, unlike the other NS parameters, e.g. spin 
period and magnetic field, NS mass is only measured in the 
binary system (e.g. Freire 2008; Lorimer 2008; Lyne & Smith 
2005). Therefore, the statistics of the measured NS masses may 
provide information about the NS accretion history in the bi- 
nary phases (e.g. Stairs 2004; Manchester 2004; Bhattacharya 
& van den Heuvel 1991). 

An accurate measurement of a NS mass in a pulsar bi- 
nary system needs five relativistic post-Keplerian parameters, 
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e.g., the periastron advance, time dilation, orbital shrinking 
rate, and Shapiro delays, which can in principle be measured. 
All of these relativistic parameters place constraints on the NS 
masses, and when three are measured, an accurate determina- 
tion of NS masses becomes possible (e.g. Lorimer 2008; Freire 
2004, 2008ab, 2009). The NS masses have been measured 
precisely in double neutron-star (DNS) systems, such as the 
first discovered pulsar PSR B 19 13+ 16 (Hulse & Taylor 1975; 
Taylor & Weisberg 1982) and double pulsar system PSR J0737- 
3039 (Burgay et al. 2004; Lyne et al. 2004; Kramer & Stairs 
2008), because the eccentric orbits of both systems provide 
well-measured relativistic parameters. Unlike DNSs, masses of 
millisecond pulsars (MSPs) are not easy to determine, since 
their binary orbits are so circular (or of such low eccentricity) 
that normally no sufficient relativistic effects can be used to 
provide extra equations to solve the masses (Freire 2000; Freire 
et al. 2004). Therefore, the masses of MSP systems are often 
measured with large errors, such as PSR J0514C4002A (Freire 
et al. 2004), except in cases of high eccentricity. The observa- 
tions of relativistic parameters in pulsar binary systems have 
presented the first application of general relativity and pro- 
vided the most widely available laboratories for testing theories 
of gravitation (e.g. Hulse & Taylor 1975; Taylor & Weisberg 
1982; Weisberg & Taylor 2003; Thorsett et al. 1993; Stairs et 
al. 2002). 

Mass measurements are also possible in X-ray binaries, 
where a neutron star X-ray pulsar and an optical companion 
reside. Careful monitoring of the cyclical Doppler shifts of the 
pulse period and Doppler shifts of the spectral features of the 
optical companion can be used to determine the orbital period 
as well as the radial velocity, which provide/infer the mass 
function of the system. Both masses are known when the in- 
clination angle of an eclipsing binary system can be measured 
(e.g. van Kerkwijk et al. 1995; Jonker et al. 2003). The accu- 
racy of this method is not so high as that of measuring DNS 
mass, usually being affected by an error of about ~ 10% or 
more (see Table 1). 

Thorsett and Chakrabarty (TC99) presented the results of 
a statistical study of 19 NS binary systems, and obtained a 
Gaussian distribution of mass around 1.35 M Q , with a narrow 
deviation of 0.04 M . The sample has increased significantly 
since then. There are now about 61 NSs with measured (esti- 
mated) masses in various types of binary systems. 

In this paper, we present a statistical analysis of the masses 
of NSs in binaries using the current data set, and investigate in 
particular the pulsar recycling hypothesis. We present a com- 
pilation of all NS mass observations in Sect. 2. In Sect. 3, we 
study the relation between the NS mass and its spin period. Our 
conclusions are given in Sect. 4. 
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Fig. 1. List of 61 measured NS masses in the different types of 
NS binary systems. Their details and references can be seen in 
Table 1-3. Vertical line M=1.4 M delineates the mass mean 
value inferred from Gaussian fitting. 




Fig. 2. Histogram of 61 measured NS masses. A Gaussian fit- 
ting curve is superimposed on the histogram plot, with the mass 
mean value 1 .40 M and standard deviation 0.18 M . 



2. Statistics of pulsar masses 

2.1. NS mass distribution 

In Tables 1-3, we list all known NSs with measured and esti- 
mated masses, including their binary parameters when avail- 
able. In Table 1, we list the 13 systems consisting of X-ray NSs 



Table 2, we first list the 18 DNSs that have masses with high 
accuracies, then 16 radio pulsars with WD companions, 3 ra- 
dio pulsars with the main-sequence star companions and one 
uncertain system. In Table 3, we also list the 10 Galactic radio 
pulsars with WD companions. 

To illustrate all NS mass distributions, a histogram of NS 



with low or high mass post-main-sequence star companions. In masses is plotted in Figf2] where a fitted Gaussian distribution 
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function is shown with a mean mass of 1 .40 M and a small 
uncertainty of 0.19 M , that is slightly higher than the previ- 
ous statistical mean value of 1.35 + 0.04 M by TC99. About 
-67% (-90%) of all NSs are within the range of 1 .2M - 1 .6M 
(l.OM - 1.8M ). The NSs with masses over 1.8M repre- 
sent about ~10% of all samples. The maximum and minimum 
values of NS masses are, respectively, 2.74 ± 0.2 M (J1748- 
2021B)and0.97±0.24M G (2A 1822-371). 

It is interesting to investigate why the present NS mass av- 
erage is higher than that measured ten years ago. The data of 
NS masses by TC99 are based on the DNSs, which are gen- 
erally less than the canonical value of 1.4 M . The present 
NS mass data includes all types of binary systems with differ- 
ent evolutionary histories. In particular, there are many NSWD 
systems, which have significantly high NS masses as shown in 
Table 1-3. 

It is generally assumed that MSPs are formed from the spin- 
up of a magnetic neutron star caused by accretion in a binary 
system (e.g. Alpar et al. 1982; Bhattacharya & van den Heuvel 
1 99 1 ; van den Heuvel 2004) . If the neutron stars were born with 
the standard pulsar type fields ~ 10 12 G, it has to be assumed 
that the field decays to ~ 10 8 ~ 9 G by accretion as well. The 
MSPs are understood to be evolutionarily linked to the long- 
lived LMXBs (e.g. van den Heuvel 2004). The evidence of a 
MSP that is linked to an LMXB was found with the discovery 
of the first accretion-powered X-ray pulsar SAX J 1808.4-3658 
(spin frequency of 401 Hz, Wijnands & van der Klis 1998). 
A consequence of the re-cycling hypothesis for the origin of 
MSPs is that the mass of a MSP should be higher than that 
of non-recycled pulsar. It has long been believed that a MSP 
should possess a higher mass than the canonical value of 1 .4 
M , e.g. ~ 1.8M , because of the significant amount of ac- 
cretion (e.g. van den Heuvel & Bitzaraki 1995ab; Burderi et al. 
1999; Stella & Vietri 1999). Thus, if this relation between MSP 
mass and accretion exists, we may expect to see it in NS mass 
statistics taken over different spin period ranges. 

We first divided all NS samples into two groups, those with 
spin periods longer than and equal to or shorter than 20 ms. 
The 20 ms dividing line was taken somewhat arbitrarily as the 
period below which a pulsar would be classified as a MSP. We 
find that the mass averages of MSPs and less recycled NSs are 
1.57 ± O.35M and 1.37 + O.23M , respectively. The expected 
trend is therefore clearly seen in the data. The above trend can 
also be seen in Fig. [3] The mass systematically decreases with 
the spin period, or alternatively, spin-up is associated with an 
increase in mass of NS. 

By dividing the pulsars into three groups, the mass av- 
erages are, respectively, M=1.57 ±0.35 M (P < 20 ms), 
M=1.38+0.23 M o (20ms < P < 1000 ms), and M=1.36+0.24 
M (P > 1000 ms). Here, we note that the average mass of 
the recycled pulsar increases with the stellar spin-up. In gen- 
eral, the spin periods and magnetic fields (B) of recycled pul- 
sars are just below the spin-up line in B-P s diagram of pulsars 
(e.g. Bhattacharya & van den Heuvel 1991; Lorimer 2008), 
where the B-P s correlation is given by P s ~ B 6 ^ 1 from the 
accretion-induced magnetic evolution model for recycled pul- 
sars (Zhang & Kojima 2006), the magnetic field and accretion 
mass correlation for recycled pulsars is given approximately 
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Fig. 3. Diagram of mass versus spin period for 39 NSs. The 
horizontal line M=l M (3.2 M ) stands for the measured min- 
imum mass (theoretical maximum mass, see Rhoades & Ruffini 
1974). The vertical line at 20 ms separates the samples into two 
groups, MSP (< 20 ms) and less recycled NS (> 20 ms). It is 
found that the mass averages of two groups are, respectively, 
1.57 + 0.35 M and 1 .37 + 0.23 M . The solid curve stands for 
the relation between accretion mass and spin period of recycled 
pulsar as described in Eq.(l) and (2), M = 1.40 + 0.43(^)" 2/3 
(M ) . 

by B ~ AM -7/4 , which infers a relation as AM ~ PJ 3 ' 2 . On the 
basis of the above estimates and arguments, we propose an em- 
pirical relation between the accreting mass (AM) of recycled 
pulsar and its spin period as 



AM = M a (P/ms) 



-2/3 



(1) 



where M a is a characteristic accretion mass when a pulsar is 
spun-up to one millisecond. The mass of recycled pulsar (M) 
increases with accretion and is roughly expressed as, 



M = M + AM 



(2) 



where Mo is the mass of NS at birth while NS spin period is as 
long as those of HMXBs. 

Exploiting Eq.(l) and (2) to fit the NS mass and spin period 
data as shown in Fig. 3, we find that M = 1.40 + 0.07 M 
and M a = 0.43 + 0.23 M . Because of the broadness of the 
initial NS mass distribution and the large errors in measuring 
NS mass, the fitting COD is as low as 0.07. 

2.2. Special DNS mass spectrum 

The mass average of all eighteen DNSs in nine systems is 
1.32 + O.14M , which is systematically lower than that of the 
less recycled NS (M=1.37+0.23 M ). The mass averages of 
the nine recycled and non-recycled DNSs are, respectively, 
1.38 + O.12M and 1.25 + 0.13M G , where the mass of recycled 
NS is generally higher than that of non-recycled one, which 
may be the indication that either the accretion induces the mass 
increase for the recycled NS or the evolution of DNS progen- 
itors makes the mass of non-recycled NS low. However, we 
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Fig. 4. Mass ratio versus orbital period diagram for 9 pairs of 
DNSs, where the vertical axis Ma/Mb represents the mass ratio 
of the recycled NS to non-recycled one. 

cannot derive how much mass is accreted into these systems, 
since for two systems (J 18 11-1736 and J 15 18+4904) both NS 
pair masses have large differences with large errors, e.g., PSR 
J181 1-1736 with 1.5+° ' 2 M and 1.06+°f M Q (see Table 2). 

The mass ratios of seven DNSs are close to unity and those 
of the other two with longer orbital periods are higher than 
unity, as shown in Figj4] It is too early to draw conclusions 
about any ratio gap, separated by the orbital period at 2 days, 
since fewer DNS samples are not sufficient to infer a warranty 
statistical result. The cause of the systematically lower mass 
values of DNS systems than the typical 1.4 M remains un- 
known. We propose that the evolution of the DNS progenitors 
may influence or interact each other, which may be responsible 
for the particular mass spectrum distributions shown above. 



2.3. On AIC mechanism for MSP formation 

Although we have focussed on the standard formation model 
(recycled NS) of MSPs which involves accretion, associated 
field decay and spin up, other models are possible (e.g. Kiziltan 
& Thorsett 2009abc). These include the often discussed pos- 
sibility of the accretion induced collapse (AIC) of a white 
dwarf onto a neutron star (e.g. van den Heuvel 1994; Verbunt 
1990; Fryer et al. 1999; van Paradijs et al. 1997; Ferrario & 
Wickramasinghe 2007). In this model, a white dwarf of mass 
> 1.2 M consisting O, Ne, and Mg (e.g. Nomoto & Yamaoka 
1992) collapses onto a white dwarf because of the accretion of 
matter during the course of binary evolution, where a NS is as- 
sumed to be born as weakly magnetic and rapidly spinning as 
those observed MSPs. 

Hurley et al. (2010) presented a comparative study of the 
expected properties of binary MSPs (BMSPs) born by means 
of NS recycling and AIC. They concluded that both processes 
produce significant populations of BMSPs that could poten- 
tially be identified with BMSPs. Furthermore, prior to the de- 
tached BMSP phase at the end of binary evolution, both the NS 
recycling and AIC binary systems may have experienced sig- 



nificant phases of accretion. Nevertheless, the AIC systems are 
likely on average to have accreted less mass. 

In Fig|3] four of twenty-two BMSPs have masses of less 
than 1.35 M , which are less than Chandrasekhar mass limit 
1.44 M Q , that may be candidate AIC MSPs. Of course, for a 
NS with initial mass of 1 . 1 M , a recycled process will also 
work by accreting 0.25 M from its companion. If we assume 
the four MSPs to be the candidate AICs, then a constraint on 
the production of AIC can be derived that no more than 20% 
(~ 4/22) of BMSPs are involved in the AIC processes. 

2.4. Pulsar: neutron star or quark star ? 

From the updated measured pulsar masses, we have insuffi- 
cient information to clearly infer the nuclear matter composi- 
tions inside the central compact objects, since we require mea- 
surements of the stellar radii to determine the nuclear matter 
properties given in FigfS] a mass-radius plot of compact ob- 
ject. Theoretically, pulsars may consist of hadronic matter only 
(Menezes & Providencia 2004a), hadronic and quark matter 
(hybrid stars) either bearing or not a mixed phase (Menezes 
& Providencia 2004b; Panda, Menezes & Providencia 2004; 
Tatsumi, Yasuhira & Voskresensky 2003) or quark matter only 
(Menezes, Providencia & Melrose 2006a; Ivanov et al. 2005). 
All calculations depend on choosing of appropriate equations 
of state based on nuclear physics and thermodynamics re- 
quirements, which enter as input to the Tolman-Oppenheimer- 
Volkoff equations. The output are a family of stars, for in- 
stance, with certain gravitational and baryonic masses, radii, 
and central energy. The maximum gravitational mass and the 
associated radius are important constraints on the equations of 
state. Generally speaking, the hadronic matter equation of state 
(EOS) produces maximum masses higher than hybrid stars, 
which in turn, give slightly higher masses than quark stars. 
Radii are usually smaller for quark stars. However, these re- 
sults are very model dependent as can easily be seen from the 
references mentioned above. 

Therefore, based on the present results we cannot deter- 
mine reliably whether the pulsar is a NS or a quark star (QS) in 
this paper. However, we note that the usage of the terminology 
NS to denote the central object of a pulsar is traditional and 
does not imply any detail of its nuclear matter composition. 

Theoretically, the NS maximum mass limit of 3.2 M was 
proposed by Rhoades & Ruffini (1974). The measured pulsar 
masses are then far below this limit, which would exclude many 
known EOS models for the behavior of matter at supra-nuclear 
densities. The possible existence of high mass NS observations 
favors a stiff EOS (e.g. Ozel 2006; on the NS stiffness see 
Stergioulas 2003). The "soft" EOS models predict lower pres- 
sures for a given density, corresponding to a less massive star, 
e.g. < 1.5M . Recycled NSs in binary systems should find that 
the stiffness increases, and that the phase transition of nuclear 
matter may occur (e.g. Glendenning & Weber 2001; Menezes 
et al. 2006b). 

The fraction of NSs with masses outside range 1.2 M - 1.8 
M is less than 20%, which would provide useful information 
about their progenitor properties in most cases. 
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System 


M(M Q ) 


M c (M e ) 


Porb(d) 


P spi „ (ms) 


eccentricity 


type 


Refs 


4U 1538-52 


1.06+°4i 


16.4 + 5J 


3.73 


5.28 x 10 5 


0.08 


HMXB 


XI 


SMCX-1 


1.05±0.09 


15.5+1.5 


3.89 


708 


< 4 x 10~ 5 


HMXB 


X2 


Cen X-3 


1.24±0.24 


19.7±4.3 


2.09 


4814 


< 8 x 10~ 4 


HMXB 


X3 


LMC X-4 


1.31+0.14 


15.6+1.8 


1.41 


1.35 x 10 4 


< 0.01 


HMXB 


X2 


Vela X-l 


1.88+0.13 


23.1±0.2 


8.96 


2.83 x 10 5 


0.09 


HMXB 


X4 




1.86+0.16 


23.8±0.2 


8.96 


2.83 x 10 5 


0.09 


HMXB 


X4 


4C/1700 - 37* 


2.44±0.27 


58+11 


3.41 


No 


0.2 


HMXB 


X5 


Her X-l 


1.5+0.3 


2.3±0.3 


1.70 


1240 


< 3 x 10~ 4 


XB 


X6 


4U 1820-30 


1 9Q+0.19 
1 '^ 7 -0.07 


< 0.106 


0.08 


6.9 x 10 5 


No 


XB 


X7 


2 A 1822-371 


0.97±0.24 


0.33±0.05 


0.23 


590 


< 0.03 


LMXB 


X8 


XTEJ2 123-058 


1 46+ 030 


53+ 28 


0.25 


3.9 


No 


LMXB 


X9 


Cyg X-2 


1.78 ±0.23 


0.60±0.13 


9.84 


No 


0.0 


LMXB 


X10 




1.5 ±0.3 


0.63 ±0.16 


9.84 


No 


0.0 


LMXB 


X10 


V395 CAR/2S 0921C630 


1.44 ±0.10 


0.35 ± 0.03 


9.02 


No 


No 


LMXB 


XI 1 


Sax J 1808.4-3658 


<1.4 


<0.06 


0.08 


2.49 


< 0.0005 


LMXB 


X12 


HETE J1900.1-2455 


<2.4 


< 0.085 


0.06 


2.65 


< 0.005 


LMXB 


X13 



* The compact object may be a black hole (Lattimer & Prakash 2007). LMXB — Low-mass X-ray binary, HMXB — High-mass X-ray binary. 
XI — van Kerkwijk et al. 1995 (M, M c , P orb , eccentricity); Robba et al. 2001 (P s ). X2 — van Kerkwijk et al. 1995 (P orb , eccentricity); van 
der Meer et al. 2005 (M, M c ); van der Meer et al. 2007 (P s ). X3 — van Kerkwijk et al. 1995; Ash et al. 1999 (P Brb , eccentricity); van der 
Meer et al. 2005 (M, M c ); van der Meer et al. 2007 (P s ). X4 — Quaintrell et al. 2003 (M=2.27,1.88M G , M c , P orb , eccentricity, P s ); Barziv et 
al. 2001 (M=1.86M Q ). X5— Clark et al. 2002 (M, M c )\ Hammerschlag-Hensberge et al. 2003 (P orb , eccentricity). X6— Cheng et al. 1995 
(P orb , eccentricity); Reynolds et al. 1997 (M, M c ); Martin et al. 2001 (P s ); van der Meer et al. (2007) (P s ). XI— Wang et al. 2010 (M, M c , 
eccentricity, P s ); Shaposhnikov et al. 2004 (M, P orb ); Dib et al. 2004 (P orb ). X8 — Jonker & van der Klis 2001 (P Brb , eccentricity, P s ); Jonker 
et al. 2003 (M, M c ). X9— Tomsick et al. 1999 (P s ); Tomsick et al. 2002 (M, M c , P orb , eccentricity). X10— Cowley, Crampton & Hutchings 
1979 (P s ); Orosz & Kuulkers 1999 (M, M c , P orb , eccentricity); Elebert, Callanan & Torres, et al. 2009a. XI 1— Steeghs & Jonker 1996; 
2007 (1.44M G ); Shahbaz, & Watson 2007 (1.370.13 M ). X12 — Elebert et al. 2009b; Chakrabarty & Morgan 1998; Jain, Dutta & Paul 
(P orb ). X13— Elebert et al. 2008; Kaaret, Morgan & Vanderspek et al. 2006 (P orb ). 




Radius (km) 

Fig. 5. NS mass versus radius plot. The EOS curves and straight 
lines follow the same meanings as those of Lattimer & Prakash 
(2004,2007) and Miller (2002), where SSI and AFO stand 
for EOSs of the quark matters. For most NSs with measured 
masses of 1.0-2.0 M Q , their nuclear matter compositions are 
difficult to distinguish as those of either neutrons or quarks, 
since NS radii cannot be precisely measured in general using 
present-day observations (e.g. Truemper et al. 2004). 



3. Summary and conclusions 

We have studied the statistical distributions of the updated mea- 
surements of pulsar masses in binary systems, and the follow- 
ing conclusions and implications are obtained: 

(1) For 61 reliably measured (estimated) pulsar masses, a 
mass average of M=1.46+0.3 M is obtained, which is higher 
than found (1.35 M ) in 1999 by TC99. 

(2) Our statistics indicate that the mass average of the more 
rapidly rotating MSPs (M= 1.57+0.35 M for P s <20 ms) is 
higher than that of the less recycled ones (M = 1.37+0. 23 M for 
P s > 20 ms). This implies that the NS masses increase in the 
accreting spin-up binary systems, while a MSP accreting about 
~0.2 M from its companion appears to be present. The rela- 
tion between the accretion mass (AM) of recycled pulsar and 
its spin period is proposed to be AM = O.43(M )(P/1to.s)~ 2/3 . 

(3) The statistics of 18 DNSs indicate that their mass aver- 
age M=1.32+0.14 M is systematically lower than the typical 
mass value of the less recycled PSRs, which seems to imply 
that the mass formation or evolution history of DNS should 
differ from those of the other binary systems. 

(4) Apart from the standard recycled processes for the for- 
mation of MSPs, the mechanism by AIC of accreting white 
dwarfs is investigated by the MSP mass distribution, since AIC 
needs the mass of MSP to be less than the Chandrasehkar mass 
limit 1.44 M . If the AIC explodes after accreting -0.1 M Q of 
crust, then fewer than 20 % of BMSPs are inferred to be in the 
AIC processes, which provide a quantitative constraint on the 
formation rates of AIC MSPs. 
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Table 2. Parameters of radio binary pulsars 



System 


M(M ) 


M C (M ) 


Porbid) 


P spi „(ms) 


eccentricity 


type 


Refs 


J1518+4904 


1 56 +a2 ° 

1-J -1.20 


1 05 +1 21 


8.63 


40.9 


0.249 


DNS 


Rl 


J1811-1736 


i c+0.12 
J -0 4 


L06 -of 


18.8 


104.2 


0.828 


DNS 


R2 


J1829+2456 


1 15 +61 

1 • 1 25 


i oc+0.46 


1.176 


41.0 


0.139 


DNS 


R3 


B1534+12 


1.33+0.0020 


1.35+0.0020 


0.421 


37.9 


0.274 


DNS 


R4 


B1913+16 


1.44+0.0006 


1.39+0.0006 


0.323 


59.0 


0.617 


DNS 


R5 


B2127+11C 


1.35+0.080 


1.36+0.080 


0.335 


30.5 


0.681 


DNS 


R6 


J0737-3039A(B) 


1.34+0.010 


1.25+0.010 


0.102 


22.7 (2773) 


0.088 


DNS 


R7 


J1756-2251 


1 40 +004 

1 o 06 


1 i O+0.06 
1 - 18 -0 04 


0.320 


28.5 


0.181 


DNS 


R8 


J1906+0746® 


1.25 


1.37 


0.166 


144 


0.085 


DNS 


R9 


J0437-4715 


1.58+0.18 


0.24+0.017 


5.74 


5.76 


1.9 x 10~ 5 


NSWD 


R10 


J062 1 + 1002 
J075 1 + 1807 


1 70 +0 - 5 V°- 32 1 

— U.Oj u.zy 

1.26 + 0.14 


n q 7 +0.43/+0.27n 
i n 9/1 V n is/ 

U.ZH U. 1 J J 

0.19+0.03 


8.32 
0.263 


28.9 
3.48 


0.003 

3 x 10-" 


NSWD 
NSWD 


Rll 

R12 




2. 1^9 j (corrected) 


0.19+0.03 


0.263 


3.48 


3 x 10-" 


NSWD 


R12 


J1012+5307 


1.7+1.0 


0.16 + 0.02 


0.605 


5.26 


< 10~ 6 


NSWD 


R13 




1.64 + 0.22 


0.16 + 0.02 


0.605 


5.26 


< io- 6 


NSWD 


R13 


J1045-4509 


< 1.48 


0.13 


4.08 


7.47 


< 10~ 5 


NSWD 


R14 


Jl 141-6545 


1.27+0.01 


1.02+0.01 


0.198 


394 


0.172 


NSWD 


R15 




1.3+0.02 


0.986 ± 0.02 


0.198 


394 


0.172 


NSWD 


R15 


J1713+0747 


1.53^(1.6 + 0.24) 


0.33 ± 0.04 


67.83 


4.75 


7.5 x 10~ 5 


NSWD 


R16 


B 1802-07 




0.36*$ 


2.62 


23.1 


0.212 


NSWD 


R17 


J 1804-27 18 


< 1.73 


0.2 


11.1 


9.34 


4 x 10~ 5 


NSWD 


R18 


B 1855+09 
J 1909-3744 


i so+0.10 
• 3 -0 13 

1.44+0.024 


ri 97+O.OIO 
— 014 

0.20+0.0022 


12.33 
1.53 


5.36 
2.95 


2.2 x 10~ 5 

io- 7 


NSWD 
NSWD 


R19 
R20 


J20 19+2425 


< 1.51 


0.32 - 0.35 


76.5 


3.93 


1.1 x 10~ 4 


NSWD 


R21 


B2303+46 


1.34+0.10 


1.3+0.10 


12.34 


1066 


0.658 


NSWD 


R22 


J0437-4715 


1.76+0.20 


0.25+0.018 


5.74 


5.76 


1.918 x 10~ 5 


NSWD 


R23 


J1023+0038 


1.0-3.0 


0.14-0.42 


0.198 


1.69 


< 2 x 10~ 5 


NSWD 


R24 


J1738+0333 


1.6 + 0.2 


0.2 


0.354 


5.85 


1.1 x 10~ 6 


NSWD 


R26 


J0045-7319 


1.58+0.34 


8.8+1.8 


51.17 


926 


0.808 


NSMS 


R27 


J1740-5340 


1.53+0.19 


>0.18 


1.35 


3.65 


< IO 4 


NSMS 


R28 


J1903+0327 


1.67 + 0.01 


1.05 


95.17 


2.15 


0.437 


NSMS 


R29 


J1753-2240 


~ 1.25 


~ 1.25 


13.64 


95.1 


0.304 


uncertain 


U 



DNS — double neutron star; NSWD — pulsar-white dwarf binary; NSMS — neutron star/main-sequence binary; ® The recycled NS should be 
the companion because of the strong magnetic field of PSR J1906+0746 ~ 10 12 G. Rl — Nice et al. 1995 (P orb ,P s , eccentricity); TC99 (M,M C ); 
Janssen et al. 2008 (m p < 1.17 and m c > 1.55M ). R2— Lyne et al. 2001 (P orh ,P s , eccentricity); Lorimer et al. 2008 (M,M C ); Breton, 2009 
(M,M C ). R3— Champion et al. 2004 (P orb ,P s , eccentricity); Lorimer et al. 2008 (M,M C ); Breton et al. 2009 (M,M C ). R4— Wolszczan 1991 
(P orb ,P s , eccentricity); Stairs et al. 2002 (M,M C ). R5— Hulse & Taylor 1975 (P orh ,P s , eccentricity); Weisberg & Taylor 2003 (M,M C ). Re- 
Anderson et al. 1990 (P„ rh ,P s , eccentricity); Jacoby, Cameron & Jenet et al. 2006; TC99 (M,M C ). R7 — Burgay et al. 2003 (P„ rh , P s , eccentricity); 
Lyne et al. 2004 (M,M C ). R8— Manchester et al. 2001 (P orb , P s , eccentricity); Faulkner et al. 2005 (M,M C ). R9— Lorimer & Stairs 2006 (P orb , 
P s , eccentricity); Kasian et al. 2007; Lorimer et al. 2008 (M,M C ); Breton et al. 2009 (M,M C ). R10— Johnston et al. 1993 (P orb , P,, eccentricity); 
van Straten et al. 2001 (M,M C ). Rl 1— Camilo et al. 1996 (P orb , P s , eccentricity); Splaver et al. 2002 (M,M C ). R12— Lundgren et al. 1995 (P orb , 
P s , eccentricity); Nice et al. 2004 (M,M C ); Nice et al. 2005, Nice et al. 2008 (M,M C ). R13— Nicastro et al. 1995 (P orb , P s , eccentricity); van 
Kerkwijk et al. 1996, 2005; Callanan et al. 1998; TC99 (M,M C ). R14— Bailes et al. 1994 (P orb , P s , eccentricity); TC99 (M,M C ). R15— Kaspi 
et al. 2000 (P orb , P s , eccentricity); Burgay et al. 2003 (M,M C ); Bailes et al. 2003; Bhat & Bailes, 2008 (M,M C ). R16— Foster et al. 1993 (P mb , 
P s , eccentricity); Splaver et al. 2005 (M,M C ). R17— DAmico et al. 1993 (P orb , P s , eccentricity); TC99 (M,M C ); Lorimer et al. 2008 (M,M C ); 
Breton et al. 2009 (M,M C ); Freire 2000. R18— Lorimer et al. 1996 (P orh , P s , eccentricity); TC99 (M,M C ); Breton et al. 2009 (M,M C ). R19— 
Segelstein et al. 1986 (P„ rh , P s , eccentricity); Nice, Splaver & Stairs 2003 (M,M C ). R20 — Jacoby et al. 2003 (P orb , P s , eccentricity); Jacoby 
et al. 2005 (M,M C ). R21— Nice et al. 1993 (P orb , P s , eccentricity); Nice et al. 2001 (M,M C ). R22— Dewey et al. 1985 {P orb , P s , eccentricity); 
Kerkwijk & Kulkarni 1999 (M,M C ). R23— Johnston et al. 1993 (P orb , P s , eccentricity); van Beveren et al. 2008 (M,M C ). R24— Archibald et al. 
2009 (P^/^eccentricit^M.M,). R26— Jacoby, PhD thesis, (2004); Freire PhD thesis, 2000. R27— Bell & Bessell et al. 1995 (M,M C ); Kaspi, 
Bailes & Manchester et al. 1996 (P orh , P s , eccentricity). R28— Kaluzny et al. 2003 (P orb ,M); D'Amico et al. 2001 (P s ,eccentricity,M c ). R29— 
TC99 {P orb , P s , eccentricity, M c ). R29 — Champion et al. 2008 (P orb , P s , eccentricity, M, M c ); Freire et al. 2009. U — uncertain companion type, 
Keith et al. 2009ab (P orh , P s , eccentricity, M, M c ). 

(5) The nuclear matter compositions of the less massive 
DNSs and heavier MSPs may be different. During accretion, 
the matter phase transition from the 'soft' EOS to 'stiff' EOS, 
or even the matter transition between the neutron and quark 



may be possible (Menezes 2006b), which would provide clas- 
sifications of the nuclear matter inside DNSs and MSPs. 

Moreover, the newly measured mass 1.97 ± 0.04 M Q of a 
MSP PSR J1614C2230 with a spin period of 3.15 milliseconds 
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System 


M(M ) 


M C (M ) 


Porb(d) 


/>,,„„ (ms) 


eccentricity 


type 


Refs 


J0024-7204I(B002 1-721) 


1.44 


0.15 


0.23 


3.49 


6.3 x 10- 5 


GC 


Gl 


J0024-7204H(B002 1 -72H) 


1 41+0.04 


n i O+0.086 


2.380 


3.21 


0.071 


GC 


Gl 


J1518+0204B(B1516+02B) 


2.08±0.19 


> 0.13 


6.860 


7.95 


0.14 


GC 


G2 


J1911-5958A 


1 40 +016 


0.18 


0.837 


3.48 


< 10 5 


GC 


G2 


J1802-2124 


1.21+0.1 


> 0.81 


0.699 


12.65 


3.2 x 10~ 6 


GC 


G3 


J1824-2452C 


< 1.367 


> 0.26 


8.078 


4.158 


0.847 


GC 


G4 


J0514-4002A 


< 1.52 


>0.96 


18.79 


4.99 


0.888 


GC 


G5 


J1748-2021B 


2.74 ± 0.2 


> 0.11 


20.55 


16.76 


0.57 


GC 


G6 


J 1 748-202 llfTer 5 I) 


1.3 + 0.02 


0.24 


1.328 


9.57 


0.428 


GC 


G7 


J1748-2021J(Ter 5 J) 


1 88+0-02 


0.38 


1.102 


80.34 


0.35 


GC 


G7 



GC — Globular cluster pulsars. Gl — Manchester et al. 1991; Freire et al. 2003; Lorimer et al. 2008. G2 — Wolszczanet al. 1989; Bassa et 
al. 2006; Cocozza et al. 2006; Freire et al. 2007; Lorimer et al. 2008; Freire et al. 2008b. G3— Lorimer et al. 2008; Lorimer et al. 2008; 
Faulkner et al. 2004; Ferdman et al. 2010 (1.24+0.1 1M Q ). G4— Ransom and Freire, 2009. G5— Freire & Ransom 2007. G6— Freire & 
Ransom, 2008; Freire et al. 2008a; Freire 2009. G7— Ransom et al. 2005. 



seems to hint that either MSP accretes more mass from its com- 
panion or a high mass of pulsar is brought in born (Demorest 
et al. 2010). 
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